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ABSTRACT 

We analyze the properties of steady and time-dependent C shocks under condi- 
tions prevailing in giant molecular clouds. For steady C shocks, we show that ioniza- 
tion equilibrium holds and use numerical integrations to obtain a fitting formula for 
the shock thickness mediated by ambipolar diffusion, Lshock oc no~^/^vo^^^-Bo^^^XiO~^- 
Our formula also agrees with an analytic estimate based on ion-neutral momentum ex- 
change. Using time-dependent numerical simulations, we show that C shocks have a 
transient stage when the neutrals are compressed much more strongly than the mag- 
netic field. The transient stage has a duration set by the neutral-ion collision time, 
iAD ~ -^shock/^^drift ~ 0.1 — 1 Myr. This transient creates a strong enhancement in the 
mass-to-magnetic flux ratio. Under favorable conditions, supercritical prestellar cores 
may form and collapse promptly as a result of magnetic flux loss during the transient 
stage of C shocks. 

Subject headings: ISM: magnetic fields — MHD — shock waves — stars: formation 



1. Introduction 

Within giant molecular clouds (GMCs), dense grayitationally bound cores form and collapse to 
create protostars (jShu et al.l 119871 : iMcKee &: Ostrikeii 120071 : 1 Andre et al.ll2009l ). Supersonic turbu- 
lence is believed to strongly affect the core f ormation and evolution p rocesses, with post-shock dense 
regions the most susceptible to collapse (see iGong Sz Qstrikeri (120 111 ) and references therein) . These 
processes can be modified significantly by the interstellar magnetic field. Sufficiently strong mag - 
netic fields, if they are well-coupled to the gas, can entirely prevent collapse (jMestel fc Spitzerlll956l ): 
this can be expressed in terms of a maxi mum ratio of mass to magrietic flux, M/$, or of surfac e 
density to magnetic fleld strength, T,/B fiMouschovias fc Spitzerlll976l : iNakano fc Nakamuralll978l ). 
However, in a partially-ionized medium, magnetic fields are coupled to the neutrals only through 
ion-neutral c ollisions. This amb ipolar drift modifies the dynamical effect of mag netic fields on the 
neutr al gas (jMouschoviasI Il979l ) , in particular altering the character of shocks (jPraine &: McKee 
1993). 
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In ideal MHD, the fluid and magnetic fields are perfectly coupled by assumption. When 
flow velocities exceed the relevant signal propagation speeds for a magnetized medium, discontinu- 
ities representing shock fronts (jump shock or J-type shock) can form. The cornpress ion ratio is 
parametrized by the particle density, inflow velocity, and magnetic field (e.g.. IShulll992l ). However, 
in lightly ionized clouds, velocity differences that would produce a J shock in ideal MHD are small 
compared to the magnetic signal speed ("Alfven speed") in the ionized medium, VA,i = B j ^f^TTpi. 
Ions and magnetic fields therefore smoothy transition between upstream and downstream conditions 
without discontinuities. As a result of the ion-neutral drag forces, the transition in the neutrals is 
also modified a nd all physica l quantities vary smoothly in the shock region, forming a continuous 
(C-type) shock (jDraindll98Cl ). In a steady C shock, upstream and downstream values of the neutral 
density, velocity, and magnetic field are the same as for a J shock. Thus, upstream and downstream 
values of the mass-to-magnetic flux ratio (per unit length parallel to the shock) are the same. Many 
studies of C shocks h ave investigated thei r formation (e . g. ISmith &: Mac Lowlll997l ). structure (e.g. 



Mac Low et al.lll995l ). and s tability (e.g. 



Wardle 



and emission properties (e.g. iDraine et al 



1990 



Stone 



19831 : IPineau des Forets et al.lll997l l. 



19971). as well as detailed chemical 



Ambipolar diffusion may play a key role in the star-forming process. In the traditional picture, 
quasi-static prestellar cores for m by gravitational ly-driven ambipolar diffusion in magnetically- 
supported clouds (see review bv I Andre et al.l (j2009l )). For a star to form out of gas that is initially 
stron gly magnetized, dense cores must lose in agnetic support so that gravitational collapse can take 
place ( Mouschovias 1978 ; Lizano &: ShulE989 ). If the magnetic pressure in a gravitationally-confined 
core exceeds that in its surroundings, the gradient in magnetic pressure makes the magnetic field 
(and ions) tend to exp and. The neut rals will be left behi nd as a supercritical core as the magnetic 
field diffuses outward (INakand Il979l ) . More realistically, iMouschoviad (119791 ) argued that a cloud 
does not need to lose magnetic flux as a whole to collapse. Rather, ambipolar diffusion redistributes 
the mass within dense clumps, with the neutrals diffusing inward while the magnetic field threading 
the outer region is left behind. The duration of the ambipolar diffusion process can be considerably 
longer (up to a factor of 10) than the gravitational fre e-fall timescale although the evolution is 
more rapid if cores are initially closer to critical (e.g. 
200lh . 



Mouschovias h Cioleklll999l : ICiolek fc Basu 



Observationally, the prestellar core lifetime can be estimated by calculating the ratio of the 
number of cores with embedded young stellar objects (YSOs) to the number of preste llar cores, 

which should be comparable to the ratio of protostar lifetime to the prestellar core lifeti me (ILee &: Mvers 

1999 ) . Several studies ha ve suggested a prestellar core lifetime of ~ 10^ yr, or (2—5)% (jWard-Thompson et al 
20071 : Evans et al.ll2009l ). This value is much lower than expected from the magnetic-dominated 
model. In addition, in the turbulence-controlled regim e where the magnetic field and ambipolar 
diffusion play minor roles (IMac Low &: KlessenI |2004|). ideal M HD simulations have shown that 
cores only live for (1 — 2)fff (e.g. IVazquez-Semadeni et al.ll2005l ). after which they either collapse 
or re-expand. This would not permit an extended period of ambipolar diffusion. 



Several studies have suggested that turbulence in GMCs can accelerate ambipolar diffusion 
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and star formation, by introducing large local gradient s and n onlinearities. Considering small-scale 
fluctuations in a background field, iFatuzzo &: Adama (j2002l ) analytically showed that turbulence 
can enhance t he am bipolar diffusion rate by a factor of 2 — 3 for typical conditions in GMCs. 
Heitsch et alj (|2004l ) investigated this problem numerically in a 2.5-dimensional geometry and con- 
cluded that th e enhanced diffu s ion ra te m ust be balanced a gainst large-scale compressive flows. 
Independently, iLi Nakamural (120041 ) and iNakamura &: Lil (120051 ) noted that the failure of the 
standard theory to predict core formation timescales indicates that dense clumps may not have 
formed quasi-statically through ambipolar diffusion. By performing two-dimensional simulations 
of magnetized sheetlike clouds, they found that with sufficiently strong turbulence, dense filaments 
can form from magnetic-field-dominated clouds in one turbulence crossing time (i ~ 10*^ yr). 



T urbulence-accelerated, magnetically-regulated star formation was studied by iKudoh &: Basu 
(|2008l ) using three-dimensional simulations, including self-gravity and ado pting hydrostat ic equi - 



librium in the vertical direction as an initial condition. More recently, iKudoh &: Basul (|201ll ) 
conducted a parameter study of fragmentation in magnetically subcritical clouds regulated by am- 
bipolar diffusion and nonlinear turbulent fiows. They concluded that the core formation time is 
strongly affected by the turbulence speed and the density in compressed re gion. These and other 
recent simulations with both strong turbulence and ambipolar diffusion (e.g.. INakamura Lill2008l ) 
are consistent with observations in terms of the cor e evolution time, the relatively low efficiency 
of star fo rmation (~ 3 — 6%, see Evans et al.l (j2009l )). and the core structure (subsonic infall mo- 



tions, see lLee et al.l (|l999l )). However, the fundamental physical process driving core formation via 
turbulence-enhanced ambipolar diffusion, as well as its dependence on environmental parameters, 
still remain unclear. 

To investigate this problem, we consider the simplest possible time-dependent problem with 
large spatial gradients: a one-dimensional high-speed converging fiow that shocks. In order to 
clearly distinguish the effect of ambipolar diffusion from other dynamics, we neglect the self-gravity 
of the gas. We also focus on the simplified case in which the infiow velocity is perpendicular to the 
magnetic field lines; more general geometry (i.e. oblique shocks) is discussed in Appendix lAl 

When gas is compressed by converging fiow, the neutrals are pushed to accumulate downstream. 
The ion density and magnetic field strength, however, will be only moderately enhanced since the 
magnetic pressure resists strong compression. These lagging ions exert a drag force on neutrals, 
reducing the streaming of neutrals into the post-shock region. The momentum exchange between 
neutrals and ions speeds up ions, increases the compression of the magnetic field, and reduces 
the post-shock density of the neutrals. Over time, a steady C shock develops. However, at early 
stages, for an interval com parable to the neutral-ion collision time, the neutrals do not exper ience 
drag forces from the ions (IRoberge &: Ciolekl l2007l : Ivan Loo et al.l l2009l : lAshmore et al.l l20ld ) . As 
a consequence, the initial shock for the neutrals is essentially unmagnetized, and the neutrals can 
be very strongly compressed. If the gravitational collapse timescale is sufficiently short, and a 
dense enough layer of gas builds up, the magnetically-supercritical region may be able to collapse 
gravitationally before a steady C shock structure forms. The transient ambipolar diffusion process 
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in shocks may help to explain the physics of turbulence-accelerated, magnetically-regulated star 
formation. 

In this paper, we first revisit the steady-state structure of C-type shocks in conditions appro- 
priate for GMCs, in particular allowing for varying ionization fraction. By fitting the results of 
steady one-dimensional solutions, we obtain an expression for the C shock thickness as a function 
of the upstream density, the velocity, the magnetic field, and the ionization fraction. These C 
shock thicknesses are comparable to, or exceed, the size of observed cores. We then consider time- 
dependent shocks, which we follow by implementing ambipolar diffusion in the MHD code, Athena. 
Our simulations suggest that under some circumstances, transient C shocks make it possible for 
a magnetically subcritical cloud to form supercritical dense cores, which would then be able to 
collapse promptly. We show more generally that the mass-to-flux ratio is signiflcantly increased by 
ambipolar diffusion in transient post-shock regions, compared to the value that would hold under 
ideal MHD or in a steady C shock. 

This paper is organized as follows. The model and the governing equations are described in 
Section [2j In Section [3] we investigate the structure of steady C shocks, and obtain (analytically and 
numerically) an explicit formula for the dependence of shock thickness on environmental parameters. 
In Section |4] the time-dependent numerical method is described, and we show that in the transient 
early development of C shocks, the post-shock ratio of density to magnetic field is very large. 
In Section \5\ we discuss mass-to-flux ratios of shocked gas, which we use to quantify the effect of 
ambipolar diffusion. A parameter study of the duration and effect of transient C shocks is presented 
in Section [6l We summarize our conclusions in Section [71 



2. Dynamical Equations and Model Parameters 

2.1. Basic Equations 

For a partially ionized medium with a drag force between ions and neutrals, the neutral 
fluid equations are 



Pn 



dpn 
dt 

dt 



+ V • (p„v„) = 0, 



+ (Vn • V) V„ 



(la) 

(lb) 



which represent conservation laws of mass and momentum, respectively. The corresponding mo- 
mentum equation for the ionized fluid and magnetic induction equation are 



Pi 



|^ + (v..V)v„ 



+ VP, - — (V X B) X B = -fd, 
47r 



c?B 

— + V X (B X V,) = 0. 



(2a) 
(2b) 
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We discuss the ion density evolution below; this must take into account ionization and recombina- 
tion. 

The ion-neutral drag force per unit volume is 

fd = apnPi (vj - v„) , (3) 

where |vj — v„| is the slip speed, and a = (crvrci) / {lJ'n + fJ-i) is the collision coefficient with the 
collisional cross-section a. The mean neutral and ion molecular weight Hn and /Xj are applied here 
so the number density is = pn / Pn, ni = Pi/ Pi. For simphcity, we shall assume an isothermal 
equation of state, Pn = c^nPm Pi = c^iPi, and = P/p = kT/p. 



2.2. Steady State One-dimensional Shock Equations 

We now consider one-dimensional solutions that are steady, d/dt = 0, in the shock frame. 
We assume the magnetic field is parallel to the shock front. The x coordinate is taken to be 
perpendicular to B and the shock front. We define the compression ratio of neutral density induced 
by the shock: 

Pn = Pnfirn, (4) 



where r„ — )• 1 upstream, and r„ — )■ const, downstream. Since PnVn = const, from Equation (|lap . 
Vn = Vnfi/rn, where Vnfi is the neutral speed far upstream. 

Since magnetic flux is conserved, ViB = const, in the gas. Far upstream, B ^ Bq = 
const, Vi — )• Vifi = const.. We define the compression ratio for magnetic field such that 

B = tbBo, (5) 

and Vi = Vifi/rs with — )• 1 upstream and rs — )• const, downstream. 

For regions far from the shock there is no structure in the fluid, dv/dx — )■ 0, dp/dx — t- 0, 
dB/dx — >• 0. For Equations (jlbp and (j2ap . this means Vi = Vn far upstream and downstream. 
Therefore Vnfi = Ui^ = vq far upstream, and Vn = vb ^ far downstream. The velocities 
of neutrals and ions are therefore given in terms of the upstream shock-frame speed vq and the 
compression ratios at any x as 

Vn = — (6) 

rn 

and 

Vi = —. (7) 
To simplify the equations, we define an ion compression ratio 



Pi = Pi^on, 



(8) 
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where pi^ is the upstream ion density, and rj — t- 1 upstream, rj — )• const, (not necessarily equal to 
rj) downstream, similar to r„ and r^. 

The momentum equations can now be expressed in dimensionless form as 



TT - + (^") = ^^"^nU - - - , (9a) 
dx \rnj ox vq \rB r„ 



-■^ ( -) + --^ (^^) + H) = -^-^^^n^^ (--- 

po rs ox \rB J Po P'i ox p Ox ^ ' vq \rB r„ 

in which M and /3 are two dimensionless parameters defined as 

\ 2 1 d2 1 / \ 2 



(9b) 



that is, upstream values of the square of Mach number and (half of) the square of the Alfven Mach 
number of neutrals, respectively. In Equations ([9]) — (jlOp and subsequently, we use the shorthand 
notation Csn — ^ Cg, Pn,o Po = Pnn-o, and VAn,o va,o- The drag force terms on the right-hand 
sides of Equations ([9a|) and (|9bp have equal magnitudes and opposite signs. Note that although 
Equations (|9ap and ()9bp represent the case with magnetic field parallel to the shock front, the 
results for the case with more general geometry are qualitatively similar (see Appendix [X] for 
detailed discussion). 



2.3. Governing Ordinary Differential Equation 

Typically, we have Pi/pn ~ 30/2.3 13, and 

(avrp]) 2 X 10^^ cm^s~^ ■? i i , ^ 

a = ^ « = 3.7 X 10^^ cm=^s-^g-^ (11 

Pi + pn 32.3m,H 



(jPraine et al.l Il983l ). The Mach number A4 is generally at least ~ 10, the plasma parameter 
is uncertain, but presumably /3 ~ 0.01 — 1, and since we are considering lightly ionized fluid, 
Xifi = riifi/no is a very small number, ~ 10"^ (here, no = po/ (2.3mH))- The compression ratios r„, 
rs, and rj are dimensionless and are maximal downstream, with typical values ~ 10. Therefore, 
the last term on the left-hand side in Equation (I9bp dominates over the other two terms. 



Retaining only the largest terms in the ion momentum equation yields 



/3^A^V„r. f J- - 1 ) . (12) 



dx vq \rB 

Using this result, the neutral momentum equation can be written as 



dx \ Vn J dx 13 dx 



or 



+ rn + -^ = const. = M^ + 1 + 



(14) 



an expression of conservation of momentum of the magnetized medium. On the right-hand side of 
Equation (HI]), we have used r„ = 1 = upstream. Equations (fT^ and represent the "strong 
couphng" approximation, in which the full magnetic force on the ions is conveyed to the neutrals, 
i.e., 

(V X B) X B 



frf = apiPn (vi - 



47r 



(15) 



(IShulll99l Equation (27.8)). 

We can solve Equation to obtain 



1 + /3 (r, - 1) ( — - 1 



1/2 



(16) 



once r„(x) is known, this gives r^(x). The compression ratio rj for both neutrals and magnetic 
field lines in the post-shock region is obtained by setting = rj = rn in Equation (jlGD . yielding 

2/3X2 



1 + /3+ (1 + /3) V 4/3A^2 



1/2- 



Note that if PM.'^ ^ 1, for a strong shock. 



^0 



In dimensional form, this is 



rj « 9.8 



_^20_\ 1/2 ^ tiQ 

lOOcm-3/ 



km/s y \/iG 



(17) 



(18) 



(19) 



Note that for oblique C shocks, rs 7^ in the post-shock region, and they both depend 
on the angle 9 between B and v. Appendix [Al provides expressions for generalized rf{0) with 
Equation ()A16P and VBjiO) with Equation (|A15p . 



Combining Equations ()16p and (I12p . we obtain an ODE for rn. The governing equation is 



dx 



where 



D 



1 + /? (r„ - 1) 



XinnoVQ. 



(20) 



(21) 



If we use c, = 0.2 km/s (T/IOK)^/^ 
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2.4. Ionization Fraction 



To solve the ODE in Equation (j20p . we need a relation between r„ and rj. In the dense 
interstellar medium, the main source of neutral ionization is cosmic rays, while ions may recombine 
in the gas phase, or on dust grains. The evolution of ion number density can be written as 



drii 



a 



gram^j^n • 



(23) 



Comparin g the orders of magnitude of t he three coefficients, CcR ~ 1 ~ 10 s ^ for cosmic ray 



a. 



gram 



10- 



IQ-^ cm - ^s"^ (lTielensll2005l . Table 4.11), and 



ionization ( Shulll992 : Draine et al. 1983 ) , Qg 

~ 10~^^ cm^s~^ when T ~ 10 K (IWeingartner Sz Drainell200ll ). In moderate-density clouds 



10 



^ — 10 ^ and Un ~ 10^ — 10^ cm ^, so we can drop the grain surface recombination 



ni/rir. 

term. The ion balance equation becomes 

9n, 



dt 



+ V • (njVj) w Ccr'^t; 



(24) 



2.4- 1- Recombination- Ionization Equilibrium 



In solving Equation (j24p . one possible approximation is to assume ionization-recombination 
equilibrium everywhere. In this case, CcR^-n ~ Ogas"-? j so that 



rii 



On 



for 



Xio = 10^ X 



CCR 



an 



(25) 



(26) 



where the coefficient x-io ~ 1 — 20 (IMcKee et al.ll2010l ). 

If we adopt Equation ()25p . then = rn , and the governing ODE becomes 

drv, ^ -Dr^J^ (—- — 

dx 1 — V fn TB 

where is given in terms of r„ by Equation 



(27) 



2.4-2. Frozen-in Magnetic Field 



Another approach to Equation (|24p is the so-called frozen-in condition (e.g., IWardld Il99d ). 
which has been applied widely. In this approximation, ionizations and recombinations are neglected, 
so that for a steady flow, UiVi = const., which implies Vi = rs- This corresponds to a "frozen-in 
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Fig. 1. — Comparison of C shock solution with different approaches to ionization. Adopted param- 
eters are no = 500 cm~^, vq = h km/s, Bq = 10 //G, and Xio = 10- "Prozen-in(l)" means upstream 
ionization is in equilibrium, and "Frozen-in(2)" means downstream ionization is in equilibrium. 
Evidently, recombination-ionization equilibrium (open circles) is an excellent approximation to the 
exact solution (solid curve). 
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field" : the compression ratio of the magnetic field is the same as the ion flow. The governing ODE 
then becomes (using Equation p^ ) 




One thing worth noting here is that in the frozen-in approximation, the ionization fraction in 
the post-shock region will be the same as in the upstream region. Since Equation (I25p must hold 
far upstream and far downstream, we must choose whether to set Xi^ based on no or rjUQ. 



2.4-3. Explicit Solution 



We can also retain all terms in the ionization-recombination equation in our numerical inte- 
gration. Using Hi = nornXi, Equation (p^ in steady state, for one dimension, yields 



dXi CCR ttgas^-O 2 , d 

-r- = XiVnTB + Xi — 

ax vq vq ax 



(29) 



Here, is given in terms of r„ by Equation (|T6|) . By integrating Equations (j29j) and (pO|) together, 
we can calculate the explicit solution for the steady C shock system. 



2.4.4- Comparison of Ionization Treatments 

To compare ionization-recombination equilibrium and the explicit solution, we choose just an 
upstream value Xio- For the frozen-in field case, we must also choose whether our solution will have 
the same upstream ionization fraction as the equilibrium case, or the same downstream value as 
the equilibrium case. Therefore there are four different cases for us to compare. 

An example comparing the shock solutions for the four different ionization choices is shown in 
Fig. [TJ Evidently, the approximation of ionization-recombination equilibrium yields a solution very 
close to the explicit solution. We have found that this is true for the full range of parameters of 
interest, no ~ 10^ to 10^ cm~'^, ~ 1 to 10 km/s, So ~ 1 to 15 /iG, XiO ~ 1 to 10. Henceforth, 
we shall adopt ionization-recombination equilibrium and use Ui oc n^^ so that = r^^, and 
Equation (p7|) governs steady C shocks. 



3. Steady C Shock Thickness 



For any given set of parameters no, vq, Bq, and XiOj Equation ([27|) can be integrated to obtain 
a steady C shock solution. However, it is also useful to obtain estimates of the dependence of 



-li- 



the C shock thickness on the basic flow parameters. This parameterization is potentially useful in 
diagnosing magnetic field strengths from observations. In addition, it provides a helpful guide to 
assessing the scales at which ambipolar diffusion becomes important in GMCs dominated by strong 
turbulence. If, by appropriate simplifications we can integrate the governing ODE of Equation (|27p 
analytically, we can obtain an approximate expression for the shock thickness as a function of no, 
vq, Bq, and Xio- Note that, since the governing equations for oblique shocks are qualitatively similar 
to the simplified case applied here, the oblique shock thickness can be approached using the same 
methods discussed in this section (see Appendix 



3.1. Exact Solution 

From numerical integrations of Equation ()27p with a range of parameters, we have found that 
T~n/fB drops very quickly at the beginning, becomes flat in the central region, then increases rapidly 
near the other edge of the shock (see bottom panels of Fig. [2] and [3]). This behavior can be used 
to define the thickness of C-type shocks. Since the minimum of r^/r^ depends on the parameters 
(see Equation ()39p below), we should ensure that our thickness definition is insensitive to this 
value. Based on these considerations, we adopt the following definition of shock thickness for exact 
numerical solutions: 



Xa 



r„/rB=0.95' 



, Xf>Xs\ => shock thickness Lgxact = — a;s • (30) 

r„/rs=0.95 \ J \ ^ ' 



Note that for some weak shocks, Vn/rs is always larger than 0.95. Therefore this definition also 
provides limitations in the parameter space to exclude shocks which are not strong and thus do not 
satisfy our strong shock analysis. 

We have integrated the shock ODE for a range of parameters, and computed the shock thickness 
according to the definition in Equation (j30p . This is the dataset of exact solutions of C shock 
thickness over a parameter grid with 10 values of no equally spaced between 10^ and 10^ cm~^, 14 
values of equally spaced between 2 and 15 km/s, 14 values of Bq equally spaced between 2 and 
15 /.iG, and 11 values of XiO equally spaced between 1 and 21. The range of C shock thickness is 0.1 
to 20 pc in this parameter range. Note that all parts of this parameter space are not necessarily 
astronomically realistic. For example, low no and high vq is unlikely to have low Xio^ so very large 
C shock thickness is not likely to be found. 

Also note that even for C shock thickness ~ 1 pc, in a real molecular cloud all the parameters 
are likely to vary within this length scale, instead of staying constant as in our models. However, 
our solutions still provide a useful guide to approximate shock thicknesses for parameters within a 
given range. 
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Fig. 2. — Exact C shock solution (solid) compared to the "zeroeth-order" estimate of Equation ([32]) 
(circles) and an improved approximation given by Equation (I40p (triangles), for parameters no = 
500 cm~^, I'D = 5 km/s, i?o = 5 //G, and Xio = 10. 
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Fig. 3. — Same as Fig. [21 for no = 500 cm ^, = 5 km/s, Bq = 10 /uG, and Xio = 10. 
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3.2. Zeroeth-order Approximation 

We consider the relative magnitudes of the terms in Equations (jl6p and (j27p . First, since 
Cs ~ 0.2 km/s whereas vq>1 km/s, in general M.'^ is a very large number, and typically M."^ ^ 
(see Equation (|18p ). Also, from Fig. [H the ratio rn/rs is small in much of the shock region. If we 
let Vn/rB <C 1, a "zeroeth-order" approximation to Equation (f27|l is 

dx ^ M"^ 
This can be integrated analytically to yield 



(31) 



3 Dx 
2M 



-2/3 



^n{x) = ( 1 - ) ■ (32) 



This "zeroeth-order" approximation to the shock structure using Equation (j32p is shown in Fig. [2] 
and[3]for two parameter sets, in comparison to the exact solution. The zeroeth-order shock thickness 
Lzcroeth is defined as x such that Vn ^ rj in Equation ([32]) . giving 

, 2M^ ( _3/2\ 2A^2 

where the second approximation assumes a strong shock, rj ^ 1. 

Substituting Equation (j2ip for in Equation ()33p . we obtain a thickness estimate in terms of 
physical parameters 

_2vq Vq 

-t^zeroeth — TT —pr , (^04j 

or in dimensional form, 

(_^)-(_-.)(M)-. (35, 

Thus, the shock thickness increases with higher upstream velocity, and decreases with higher up- 
stream neutral density and ionization fraction. In this "zeroeth-order" approximation the shock 
thickness does not depend on the upstream magnetic field strength. From the examples shown in 
Fig. [2] and El we can see that although the zeroeth-order solution follows the general behavior of 
C shocks, it is not accurate for strongly-magnetized cases (Fig. [3|). Compared with the dataset of 
exact solutions discussed in previous section, the RMS value of (Lexact — -^zeroeth)/-^exact is 0.355, 
and the range of (Lexact - -^^zeroeth)/-^-exact is -0.8 to 0.28. 



The dependence on the velocity, ion density, and collision coefficient in Equation (I34p can be 
understood in terms of the drag force between ions and neutrals. The total momentum flux in 
neutrals entering the shock is PqVq^ . The mean drag force per volume is ~ apopifiVQ. The ratio of 
these quantities, which is the characteristic distance over which momentum exchange takes place, 
is 

L r . oc vqUq Xio • (36) 
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This dependence is similar to Equation (3.12) in lDraine &: McKed (|l993l ) if the Alfven speed in the 
fluid is similar to the upstream velocity, va ^ vq. Although they obtained an estimate using different 
assumptions and approximations, the basic idea that the momentum transfer rate determines the 
shock thickness is similar. 



3.3. Magnetic Field Influence 

To obtain a more accurate estimate of the C shock thickness, we return to the differential equa- 
tion (|9a|) for neutral momentum flux, making use of Equation (|2ip and the ionization equilibrium 
condition = r^^, 

We integrate this equation, using constant values on the right-hand-side 

,1/2 



where the minimum value of Vn/rs can be derived explicitly from Equation ()16p as 

_ 3\/3 

This yields a quadratic for function of x: 

2 



(38) 



(39) 



r 



M^ + l-D (r]/^) (l - rn + M^ = 0. (40) 



Solving Equation ()40p for r„(x) gives us another analytical approximation of the shock structure. 
When compared with the explicit solution and the zeroth-order approximation in Fig. [2] and [3l we 
can see that this correction is necessary only when the background magnetic field is strong (Fig. [3]). 

For Equation ()40p the magnetically-corrected estimate of the shock thickness {x = Lgst such 
that r„ = rf) can be written as 

i,, = (■^;- ^/)"-/-" . (41) 

D (ri ) (1 - Tn/rB) Vf 
Assuming ^ rj ^ 1 and {rn/rB)mm ^ 1, and using Equation (f2T]) . we have 

^^^^ - 77172 172- 

Note that a similar result can be obtained for the generalized case with an oblique C shock (Equa- 
tion ()A19P ). See Appendix lAl for detailed discussion. 
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Taking the strong-compression limit ~ -v/2^^o/^A,o of Equation (jlSp . we have 

^est 



97/4, 1/2 1/2 

^ ^0 ^A.O -3/4 1/2 nl/2 _l 



or in dimensional form 

Lest = 0.22 pc X 



no 



V 100cm 



-3 



-0.75 



0.5 



km/s7 V/"Gy VIO 



(43) 



(44) 



Compared with Equation (j34p . the shock thickness still depends positively on inflow velocity and 
negatively on upstream density and ionization fraction, but now a dependence on the magnetic 
field enters as well. Compared with the dataset of exact solutions discussed above, the R MS value 
oyXexact - -Lest)/L e xa,r.t if 0.13, and the range of (Lexact - Lest)/Lexact is -0.21 to 0.26. IWardle 
(|l990l ) and iLi et al.l ( 20061 ) find Lshock ~ V2^^A,o/ (api.o) in the case where ions are frozen in; this 
is smaller than Equation (f43l) by a factor 2~^/^ {vAfl/vo)^^"^- 



3.4. Numerical Approach 

Using the dataset of exact solutions discussed in Section [3TT| we construct a simultaneous linear 
fit for logLexact to log no, log^o, loguo, and logXio- We find 

n vo / \ 0.54 / D \ 0.46 , 

Over the parameter grid, the RMS value of (Lexact — -^fit)/^cxact is 0.08, and the range of (Lcxact — 
Lfit)/Lexact is -0.29 to 0.22. 

The result in Equation ()45p agrees with our expectation that the shock thickness depends on 
the magnetic field. Also, the dependences on all parameters are extremely close to Equation ()44p . 
Table [1] lists a set of model parameters (to be used in time-dependent simulations) and the C 
shock thickness based on the analytic estimate in Equation (j42p and the multivariate fit in Equa- 
tion ()45p . in comparison with the results from explicit integration of the ODE (Equation (|27p ). 
Both approaches are useful to estimate the shock thickness. 



C Shock Formation 



4.1. Numerical Algorithm for Ambipolar Diffusion 



To investigate h ow C shocks dev elop in time, we use a modified version of the numerical 
MHD code, Athena (jStone et al.ll2008l ). Athena employs a single-step, directionally unsplit Go- 
dunov scheme to obtain conservative, second-order accurate solutions of the ideal MHD equations 
([Gardiner Stonelbood l. 
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Table 1. Steady C shock Thickness Comparison 



iviocLei 


no 






-Do 


XiO 


exact 


^shock (pc) 

est. 


ni 




(cm ) 


(km/ s) 






eq. \\z i\) 


eq. [|4z|i 


cq. l|4o|l 


NOl 


100 


5 


10 


5 


3.03 


3.15 


2.89 


N03 


300 


5 


10 


5 


1.20 


1.38 


1.30 


N05 


500 


5 


10 


5 


0.82 


0.94 


0.89 


N08 


800 


5 


10 


5 


0.58 


0.66 


0.63 


NIO 


1000 


5 


10 


5 


0.50 


0.56 


0.54 


V04 


200 


4 


10 


5 


1.41 


1.68 


1.54 


V06 


200 


6 


10 


5 


1.55 


2.05 


1.92 


V08 


200 


8 


10 


5 


1.79 


2.37 


2.24 


VIO 


200 


10 


10 


5 


2.08 


2.65 


2.53 


V12 


200 


12 


10 


5 


2.38 


2.90 


2.79 


B02 


200 


5 


2 


5 


0.92 


0.84 


0.83 


B04 


200 


5 


4 


5 


1.08 


1.18 


1.14 


B06 


200 


5 


6 


5 


1.26 


1.45 


1.38 


BOS 


200 


5 


8 


5 


1.46 


1.68 


1.57 


BIO 


200 


5 


10 


5 


1.66 


1.87 


1.74 


B12 


200 


5 


12 


5 


1.89 


2.05 


1.89 


B14 


200 


5 


14 


5 


2.12 


2.22 


2.03 


XOl 


200 


5 


10 


1 


8.32 


9.37 


8.71 


X06 


200 


5 


10 


6 


1.39 


1.56 


1.45 


XIO 


200 


5 


10 


10 


0.83 


0.94 


0.87 


X15 


200 


5 


10 


15 


0.55 


0.62 


0.58 


X20 


200 


5 


10 


20 


0.42 


0.47 


0.44 
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In the strong coupling limit, the drag force = apiPn (vj — v„) is equal to the Lorentz force 
fL = [(V X B) X B] / (47r). The momentum equation for neutrals is thus identical to that in the 
ideal MHD limit. The mass conservation equation for neutrals is also the same as for ideal MHD. 
In this approximation, v, = v„ + [(V x B) x B] / {AnapiPn), so that the induction equation ()2bp 
becomes 



c?B 

— V X (Vn X B) = V X 

dt V n / 



((V X B) X B) X B 



(46) 



AirpipnO 

With v^; — Vj — the drift velocity between ions and neutrals, we can write the correction term 
in Equation (06]) in terms of a "drift" EMF, 



= Vd X B 



[(V X B) X B] X B 

Anpipna 



(47) 



In our simplified 1-D problem, B = Byy, v = v^i, and the discretized magnetic field corrected 
by ambipolar diffusion before each step, at interface position iAx and time nAt, is 



B 



in.+ l 



y\i 



in At / ,„ 



and 



y\i+i 



B„ 



B^ 

y 



Ana 



Ax 



PiPr, 



This term is implemented in Athena as an operator-split update to the magnetic field. The mesh 
resolution is set to be 0.01 pc. 



In setting the timestep, we implement the super-timestepping approach as described bv lChoi et al 



(|2009l ). choosing the factor z/ = 0.2, and taking the ambipolar diffusion timestep 

PnPi 



AtAD = 27ra (CFL number) (Ax) ■ min 

where the CFL number is set to be 0.8 in all simulations. 

For figures presenting numerical results, n„ — )• n and v„ — >• v. 



Uy 



(48) 



4.2. Convergent Flow Test 

4-2.1. Simple Convergent Flow Test 

One way to produce shocks in a numerical simulation is to use a simple convergent flow, in 
which the initial conditions are 



n = const., By = const., 



I'infiow, left half 
-■Winflow, right half. 



(49) 
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Fig. 4. — Transient C shock structure (squares) compared with ideal MHD shock (thin hnes), 
generated from convergent flow, for parameters no = 200 cm~^, vq = 1 km/s, Bq = 2 fiG, and 
Xio = 10. The central peak in density n and mass-to-flux ratio n/By is a signature of early C shock 
structure. 
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This will evolve to a dense post-shock region in the center, with outward-propagating reverse C-type 
structures at the left and right (Fig. . 

The C shock structures seen at x ~ 1.5 and x ~ 3.5 in Fig. U] are the same as the steady 
solutions obtained by integration of Equation (I27p , as confirmed by comparing the detailed profiles 
(not shown). Note that the mass-to- flux ratio n/By is analogous to r-n/rB (except not normalized 
by upstream values). The dips in n/By at the C shock locations correspond to the "well" in Tn/rs 
seen in Fig. [2] and [3l However, the central peaks in both neutral density n and the mass-to-flux 
ratio n/By in Fig. [J] are not a feature of steady C shocks. As the solution shown in Fig. [Revolves 
further in time, these peaks disappear. Thus, these peaks are a signature of transient C shock 
development, as we discuss further below. 



4-2.2. Colliding Clouds 

The initial conditions for the simple converging flow are somewhat artificial, in that only the 
velocity is discontinuous. Thus, we would like to test whether shocks formed under more realistic 
conditions also show the transient peaks in n and n/By described above. 

We consider the collision of two idealized clumps inside a large molecular cloud. We suppose 
that the two clumps are both denser than their surrounding, but the mass-to-flux ratios are the 
same throughout the whole cloud. We imagine that the large-scale turbulence in the molecular 
cloud imposes velocities such that the two dense clumps collide with each other, producing a shock. 
We simulate the scenario described, setting the background density to be 5% of the value in the 
dense clumps, and the initial velocity of this gas to be zero. We focus just on the collision region, 
so that the right and left sides of the domain are set to "clump" conditions, as in Equation (|49p . 

When the two dense clumps meet each other, a strong shock forms (Fig. [5]). Since all fluid 
variables (n. By, Vx) are smooth and continuous prior to shock formation, the features produced are 
not a consequence of discontinuous initial conditions. This test case eventually evolves to profile 
similar to that in the simple convergent flow test (Fig. [6|. The central peaks in density and mass- 
to-flux ratio show up as well. Subsequent evolution leads to a decline in the central peak in n and 
n/By (Fig. ED. 



4-2.3. Transient C shock Development 



Peaks in density above the "steady" shock solutior i have also been o bserved in other ambipolar 
diffusion simulations using different MHD codes (e.g.. IChoi et al.ll2009l ). In addition, similar tran- 
sient behavioi_of^C_slwcks_J]^s_b^ chemistry implemented 
(e.g. IChieze et al.l 119981 : Ivan Loo et al.l l2009l : lAshmore et al.M2010l ). Physically, we believe these 
peaks arise because the neutrals are effectively "unmagnetized" when the shock first forms. As 
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Fig. 5. — Two dense clumps collide with each other and produce a shock (squares). Conditions 
at a time 0.88 Myr prior to the collision are shown as thin lines for comparison. Parameters are 
no = 500 cm~^, vq = 1 km/s, Bq = 4 /xG, and Xio = 10- 
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Fig. 6. — Shock structure generated from the colhding clumps (squares), at a time 1.08 Myr after 
the stage shown in Fig. [5l Note that the central peaks in n and n/ By are qualitatively similar to 
those in Fig. [H These central peaks will then expand and smooth out (thin lines show solution 
after an additional 7.67 Myr). 
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Fig. 7. — Ptot (solid), Pb (dashed), and Pgas (dotted) in the colliding clump simulation, corre- 
sponding to the earlier (left) and later (right) stages in Fig. [6l Note that the magnetic pressure 
dominates the post-shock region, with the central peak in thermal pressure compensated for by a 
reduction in the magnetic pressure. Over time, the thermal pressure peak and magnetic pressure 
valley decline due to diffusion within the post-shock region. In the frame of the (right- or left-ward) 
expanding shock fronts (not shown), the total pressure in the post-shock region is the same as the 
total pressure upstream. Units for pressure are [2.3mH cm^^] [l km s^^] . 
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a consequence, the neutrals can be very strongly compressed, forming what is seen as a central 
density peak in Figs. HHSl 

The magnetic field, however, does not follow the initial strong compression of the neutrals. 
Instead, the overly-compressed neutrals generate higher pressure in the central regions, inhibiting 
the magnetic flux from getting in. Fig. [7] shows the total pressure Ptot = PnVn^ + -Pgas + Pb of the 
system, where Pgas = Pries'^ and Pb = B'^/Stt, with Vn measured in the "laboratory" frame. These 
three terms correspond to A1^/r„, r„, and rs^ / respectively, in Equation (fH|) . Since Athena uses 
the conservative form of the momentum equation {d{pv) /dt + dP^ot/dx = 0), Ptot must become 
constant in the post-shock region at late times. For strong shocks, the magnetic pressure term 
dominates at late times in the post-shock region. At early times, there is a slight depression of the 
magnetic field strength at the center of the shock, in order to balance the extremely high neutral 
gas pressure in the density peak. 

Combining the strong neutral compression and slight magnetic exclusion, the mass-to-flux 
ratio is elevated in the center when a shock forms. The collisions between neutrals and ions will 
gradually slow down the incoming neutrals and compress ions and magnetic field to the center. 
Meanwhile, the neutrals in the central peak diffuse outward in order to balance the increasing 
magnetic pressure and keep the total pressure constant. Eventually, the ions and neutrals interact 
sufficiently that a steady-state C shock structure develops. The post-shock n/By is the same as 
the upstream value. However, the ambipolar diffusion process takes time, and during the transient 
stage, a region of very strongly compressed neutrals will be present. 

Our finding that there is a transient stage of very strong density compression, with an enhanced 
ratio of n/By or mass-to- magnetic flux, suggests that the very early stage of shock development in 
GMCs may be particularly important to star formation. The following sections examine this idea 
further. 



5. Criticality of Clouds 
5.1. Mass-to- flux Ratio 



The mass-to-flux ratio is a crucial parameter defining whether the magnetic field can support 
a cloud against its own self-gravity. The critical value of M /^b for an uniform, spherica l cloud 
has been derived to be M/^^l^it = cq>/VG 0.126/VG (jMouschovias Spitzeil Il976l l. The 
numerical coefficie nt differs with the geom etry of the cloud: an infin ite sheet-like cloud has 
c$ = l/27r 0.16 (|Nakano k NakamuralfigTa ) . while iTomisaka et al.l dlQsd) fou nd Cds = 0.17-0.18 
for clouds with various M/^b distributions (see review bv iMcKee &: Ostrikeri (j2007l )). Since the 
value of Cj) varies o n ly ~ 10% with geometry, we choose the commonly-used c$ = l/27r (e.g. 

as a reference value, while keeping in mind 



Kudoh fc Basu 



2011 



Vazquez-Semadeni et al 



2011 



that core geometry is not explicitly defined for our slab system. 



Practically, for magnetic field in the y-direction the ratio can be written as 



M _ J pdx ■ LyLz _ Ly J pdx ^ T, 

'¥ ~ J Bydx ■ ~ jBydx jBy}' ^ ' 

where we assume that if a core formed in the post-shock region, its effective length in the y-direction, 
Ly, would be comparable to that in the x-direction, L^, so that {By) = J Bydx/ J dx ^ J Bydx/Ly. 

The mass-to-flux ratio, in units of the critical value M/^b\^j.^^ = ^27r\/G^ , is 



{By) ' Vl02icm-2 J \ pG 

To convert the column of neutrals in our simulation to A^(H), we use n = +^^116 = O.Gnn- Note 
that the true value of the normalized mass-to-flux ratio would differ from Equation (15ip by a factor 
Ly/Lx, which could be up to ~ 2. 

If the mass-to-flux ratio of a prestellar core is larger than the critical value (F > 1), i.e., the 
gravitational force exceeds the magnetic support, the core is supercritical and is eligible for collapse 
(subject to support by thermal pressure). In contrast, a subcritical core has a mass-to-flux ratio 
smaller than the critical value (F < 1), and cannot collapse unless it loses magnetic energy in 
either the strong-gravity mode (the field lines diffuse outward through ambipolar diffusion while 
gravity holds the gas material together) in which F ~ 1 is required, or the magnetic-dominated 
mode (neutral mass moves toward the center under the gravitational pull while ambipolar diffusion 
allows the magnetic field lines to remain stationary) so the mass-to-flux ratio increases. 



5.2. Bonnor-Ebert Sphere 



A typical low-mass prestellar core has Score ~ 1 Mq/ (0.1 pc) vr 
core with i?initiai ^ 27r\/G • S 

core ~ 10.7 /iG may be subcritical. 



0.007 g • cm 2, so that a 



More precisely, we consider the Bonnor-Ebert sphere radius for a core whose mean density is 
equal to the post-shock density pf, 



Rbe 



2.7c, 



(ATrGpf) 



1/2 



(52) 



(e.g. iGong &: Ostrikeij|2009l ). which is the largest sphere that can be supported by its own internal 
thermal pressure. 

We note that, from Equation (j42p . the ratio of the shock thickness to the diameter of a Bonnor- 
Ebert sphere at the post-shock density is 
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A converging flow bounded by C shocks has breadth at least twice the shock thickness. Under 
conditions in GMCs where vq/cs > 10, and Xio ^ 10, this implies that shocks are sufficiently broad 
that Bonnor-Ebert spheres can fit within the post-shock region. Thus, if magnetic fields are weak 
enough, cores could grow and collapse in post-shock gas. 



The mass-to-flux ratio for a sphere of radius i?BE in a post-shock magnetized region, without 
ambipolar diffusion, is 



M 



(54) 



with corresponding 



BE 



BE 



M/^> 



B\ 



where 



SttVG^ Pf 1.8c, 

5 -ftBE^ = 

3 Bf VAJ 



B)/ {A^pf) = rfV^Q 



1.8c, 

1/2 

r/ VA,o 



(55) 



using the shock jump conditions. Note that for a strong shock, rj « \/^vo/va,o, so that Fbe 
1.5c,/ (■uo^^A.o)^''^, 



or 



BE 



18 (-^ 



no 



V 100cm- 



0.25 



^^0 



km/s 



-0.5 



pG 



-0.5 



lOK 



-0.5 



(56) 



If Fbe is larger than 1, a post-shock region of radius ~ i?BE is dense enough to gravitationally 
collapse whether or not there is ambipolar diffusion. Otherwise, ambipolar diffusion would be 
needed for a region of size ~ iiBE to become supercritical. For the set of shock models we are 
studying (see Tabled] for inflow parameters), calculated radii and mass-to- flux ratios for Bonnor- 
Ebert spheres under post-shock conditions without ambipolar diffusion are listed in Table [2j In 
all cases, Fbe is much smaller than 1, which means no collapse at the BE scale could happen 
in the post-shock region without significant ambipolar diffusion. More generally, since Fbe ~ 
c,/ (vQVAfl)^^'^ ^ 1 under GMC conditions, most post-shock regions are sufficiently magnetized 
that gravitational collapse of low-mass cores would be prevented unless ambipolar diffusion occurs. 
Note that if p > rjpo, as is true in the candidate core material for transient C shocks, Rbe will be 
lower than the value in the table. 



6. Core Forming Process 

Current theoretical and observational work suggests that shocks produced by supersonic tur- 
bulence play a role in compressing gas to form prestellar cores. Our findings that the neutrals 
are compressed more than the magnetic field during the early stages of shock formation raise an 
interesting question: Is it possible for a subcritical cloud to form supercritical cores in shocks, which 
can then gravitationally collapse promptly? 
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6.1. Evolution of Overdense Regions 

For a prestellar core to collapse, the region must be dense enough so that self-gravity overcomes 
the magnetic support. Although self- gravity is not included in the present models, we can make an 
initial assessment of whether transient C shocks are likely to affect the ability of cores to collapse 
promptly after they form. 

In the context of our converging flow test, we shall define "candidate core" material to be 
regions where 

^ > 1-2 X ; (57) 



V -^y / background 

the background has uniform n/ By, so the time evolution of this candidate core material is easily 
calculated. Physically, this candidate core material corresponds to that in the central peak of n/By 
as shown in e.g. Fig. |4] or Fig. [HI 

For steady shocks with compression factor Vf produced by a two-sided converging flow with 
inflow speed ^inflow from both sides, the upstream speed in the shock frame is vq = ^inflow?'// {ff — 1), 
and t'shock = inflow/ i^f ~ The rate at which the column density grows for a steady shock is 
therefore 

= 2nH.inflow = 1.05 X 10 cm Myr (^00^) V W^J ' 
where n = rin^ + nne = O-Qnn is assumed. 



(58) 



For a steady shock, the growth rate of the post-shock column is independent of a, the collision 
coefficient between neutrals and ions. Fig. [8]shows evolution of the "candidate core" column density 
with different values of a, for a model with no = 100 cm^'^, ^inflow = 3.87 km/s, Bq = 10 fiG, 
and Xio = 5. In the very beginning, all of the post-shock material has n/By greater than the 
background value, because the core grows by unimpeded motion of the neutrals which do not "see" 
the ions. Thus the column of "candidate core" material initially follows Equation (I58jl . with slope 
AiV(H)/At « 4 X 10^^ cm~^Myr"^ (also shown in Fig. [8|), independent of a. 

It is evident that at some point the growth rate of the "candidate core" column decreases. 
Physically, the growth rate decreases as ions are pushed into the column by inflowing neutrals, 
which causes the magnetic flux in the "candidate core" region to increase more rapidly than the 
neutral column, and the mass-to-flux ratio to decrease. Therefore, we might expect the growth of 
the demagnetized column to slow down on a timescale ^ad ~ -^shock/'^^drift! the time for neutrals 
to travel across the shock front under the influence of ions. In timescales short compared to fAD? 
neutrals which have arrived at the center were moving fully or partly free from collisions with ions. 
These neutrals thus contribute to the column with high mass-to- flux ratio. After tAD) neutrals 
which have interacted strongly with ions dominate, and the growth rate of the low-magnetization 
column starts to decrease. 



This also corresponds to the timescale for steady C shock structure to develop, and for the 
fronts surrounding the shocked layer to expand. After this time, neutrals must then travel a greater 
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Fig. 8. — Dependence of "candidate core" material column density on the ion-neutral collisional 
coefficient, for a = ao = 3.7 x 10^^ cm^s~^g~^, a = 0.5ao, and a = 2ao. Also shown (straight 
dotted line) is the "kinematic" growth rate dN{R)/dt = 2nH'Uinflow for a steady shock. Parameters 
for this model are no = 100 cm~^, vo = 5 km/s, Bo = 10 fiG, and Xio = 5. 
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distance through the condensed magnetic field and ions to stream into the "candidate core" area. 
From Fig. [HI the "saturation" time varies, roughly inversely with the collision coefficient a. With 
Equation or (|^2|) we have Lshock (and hence tAo) oc 1/a, consistent with our numerical results. 



6.2. Time Scale and the Mass-to-flux Ratio 

As discussed earlier, our simulation results indicate that during an initial transient period, 
the neutrals are compressed much more strongly than the magnetic field. Up to a certain time, 
corresponding to the ambipolar diffusion time scale, the column density of gas with elevated n/B 
grows. After this time, the profile transitions to that of a steady C shock, with n/B equal to the 
upstream value. 

The ambipolar diffusion time scale should be comparable to the time it takes for neutrals to 
travel through the thickness of a C-type shock under the influence of ion drag. Therefore we have 

), (59) 



, -f'shock -^'shock / I 1 

IAD 



(^^drift) ^^0 \ rn TB 

where we have used I'drift = \vi — Vn\ = j^'n ^ ~ ''^"'^1 (^^e Equations ([6]) and ([7])). Assuming that 
Tn ^ fB over the shock region because a steady-state C shock has not yet formed, and using 
(tb) ~ ^//2 as an average value, we obtain 

, / \ -^'shock ^ f ^shock inn\ 

tAD « {rB) ~ (60) 

VQ 2 VQ 



In the second line, we use the estimate of Equation (|42|) for the shock thickness, and Equation p8p 
for rj in a strong shock. Note that a similar formula for more general cases with oblique shocks is 
given in Equation ()A20p . In dimensional terms, using the analytical approximation Equation (|44p 
to the C shock thickness Lshockj Equation ([60]) gives 

1„. f (Si)-"-' (m-' (62) 



VlOOcm-3/ \km/s J \lJ.G J V 10 



The time tAD can be compared to the gravitational free fall time 

-16 



/ 37r Q f n \-i/2 



to give 

+ _ / ^ \ 1/2 / „, \ 1/2 

Xio^^- (65) 



4. / \ 1/2 / \ 1/2 



ts{p) \PoJ \VA,0 
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The post-shock gas has density pj = rjpo with rj ^ -v/2 (i;o/va,o) (see Equation ([H])), which means 
that tAu/tsipf) ~ {vo/vA,o) Xio'- During the transient stage, pt > Pf, so % (pt) < tfi (p/), implying 
tAD/ts (pt) > (vo/vAfl) Xio^ ■ Thus, for strong shocks (wo/fA,o ^ 10) and low ionization conditions 
{Xio ^ 10), the transient duration ^ad will exceed the time ts (pt) for post-shock perturbations to 
develop into collapsing cores. 

If the growth rate of the neutral column density is dN^/dt ~ 2n}ivo (see Equation (j58p ). the 
final (maximum) value of the mass-to-flux ratio should be 



M 



final 



lAmudNu/dt x ^ad _ '^PoVq x ^shock?^// C^vq) _ Po-^shock ^gg^ 

^final ^f^O 



This estimate of the mass-to-flux ratio inside the pre-collapsing core depends only on the upstream 
density and magnetic field, and the steady-state C shock thickness. This can be evaluated (using 
Equation (jS])) to give Efinai = 27r\/G (M/$B)finar 



0.14 



VAflJ ^10/ 



Note that true F would differ by a factor Ly/L^ from Equations (f67|) and (j68|) : i.e. up to factor 
larger. 



Combining Equations ([61]) . (j64l) . and (j66l) . we have 



32x1/2 
67r2 



iAD = ( — ) Tfinaitfr (po) = 0.74 Tfinaitfr (po) • (69) 



Thus, shocks that are able to reach Egnai ~ 1 through transient ambipolar diffusion will do so on a 
timescale comparable to the gravitational time (po) of the large-scale cloud. Since the large-scale 
dynamical timescale (~ %(po) for a self-gravitating cloud) determines the correlation time of the 
flows that create shocks, this means that shocks will be sustained long enough for diffusion to occur. 

If ^^0 ^ ^'A.o and XiO ~ li from Equation ([68]) the candidate core will have Efinai exceeding 
unity. In this situation, a core would be able to collapse promptly, without an extended period 

1 /2 

of ambipolar diffusion, since tAD/%(pt) ^ {vo/vAfl) Tfinai > Efinai. In GMCs , the ionization 



fraction is dependent on chemical processes, with XiO ~ 1 — 20 (jMcKee et al.ll2010l ). The turbulent 



flow speed will not exceed ~ 10 km/s under realistic conditions, and no ~ lO'^ — 10^ cm~^, typicall}0. 
Therefore, Tgnai will exceed 1 (see Equation ([H7|) ) only if the upstream magnetic field in the direction 
parallel to the shock front is moderate, probably < 10 pG. The line-of-sight magneti c field stre ngths 



in molecular clouds with density < 10 cm , however, can vary over ~ 5 — 25 pG (|Crutcheiill999i : 



Keep in mind that some combinations of parameters are not astronomically realistic; e.g. high vo is unlikely to 
have low XiO, and high no is unlikely to have low Bq. 
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Crutcher et al. If the total magnetic field strength (which is always > -Blos) exceeds ~ 20/iG, 

then in order to reach Ffinai close to 1 so that pre-collapse cores can develop efficiently, converging 
flows with vq aligned < 30° to Bdoud are favored. 

In addition, we note that for itsd, cloud and f a, cloud the 3D turbulent velocity dispersion and 
mean-field Alfven speed in a cloud, a gravitationally-bound (or virialized) cloud has 

'^3D, cloud 

i cloud ~ , (70J 

VA, cloud 

SO that 

r,..,~f (^)'7^)"V;,'L. (n) 

XiQ \<^SD, cloud/ V -°0 / 

The strongest shocks will have vq ~ (730^ cloud- These regions will be able to reach Tfinai ~ 1 if the 
cloud is sufficiently supercritical (Fdoud ^ 1)) the ionization fraction is sufficiently low {xio ~ 1); 
and/or the magnetic field parallel to the shock front is weaker than the mean field threading the 
cloud (-Bcioud/^o > !)• Again, with realistic Xio and Tdoud) the most favorable circumstance for 
ambipolar diffusion to yield Tgnai > 1 is if the inflow vq is aligned locally towards Bdoud so that 

-Bcloud/-Bo >li 

Even if post-shock regions are subcritical, transient ambipolar diffusion significantly increases 
the mass-to- flux ratio compared to the value that would hold in ideal MHD. A measure of the 
importance of this effect is the ratio between Tgnai and Fbe in the post-shock region. From Equa- 
tions ([6711 and ([56]) . 

is predicted. The turbulent motions in clouds can achieve M ~ 50. With Xio ~ 1 — 20, a significant 
enhancement in the mass-to-flux ratio can be expected due to transient ambipolar diffusion. 



6.3. Simulation Results 

The estimates of Equations ()62p and ()67p can be compared to the ambipolar diffusion time 
and mass-to-flux ratio as measured directly from time-dependent numerical simulations. Examples 
showing evolution of the measured F for several different parameter values are shown in Fig. [9l 

To read the ambipolar diffusion time scale from simulations, recall that the growth rate of the 
mass-to- flux ratio inside the core decreases at time ~ i ad • We adopt a definition of t ad as the time 
when the slope of the F vs. time curve drops to 25% of its maximum value. For each simulation, we 



^We have investigated oblique shocks with nonzero B^fi ~ Bcioudcos^ = Bo cot 9 in Appendix \Kl where Bo = 
Byfi — -Bcioud sin 61 is the magnetic component parallel to the shock front. Equation ()A2ip gives an approximation of 
Tflnai as a function of 9 (the angle between Bdoud and vq). Since there is no strong dependence of Ffinai on 9, our 
1-D results (Equations (|67|l . (|68|) . and ((71}) are applicable in most cases with nonzero Bx,o- 
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Table 2. Results for Transient Mass-to- flux Enhancement 



tAD Ffinal ^corc BE SphcrG 







1 COUilr 




2tAD 




-ftBE 
(nr) 


A BE 


^ 2tAD / 
A BE 


NOl 


1.42 


1.12 


0.55 


0.37 


0.58 


0.46 


0.12 


3.08 


N03 


1.10 


1.11 


0.74 


0.64 


0.34 


0.20 


0.15 


4.27 


N05 


0.98 


1.09 


0.85 


0.78 


0.27 


0.13 


0.17 


4.59 


N08 


0.88 


1.07 


0.96 


0.95 


0.23 


0.09 


0.19 


5.00 


NIO 


0.83 


0.90 


1.02 


1.07 


0.21 


0.08 


0.20 


5.35 


V04 


1.07 


0.94 


0.59 


0.44 


0.40 


0.31 


0.16 


2.75 


V06 


1.33 


1.42 


0.73 


0.60 


0.42 


0.25 


0.13 


4.62 


V08 


1.56 


1.86 


0.86 


0.76 


0.42 


0.21 


0.11 


6.91 


VIO 


1.75 


1.88 


0.96 


0.91 


0.43 


0.19 


0.10 


9.10 


V12 


1.94 


2.13 


1.06 


1.06 


0.43 


0.17 


0.09 


11.78 


B02 


2.88 


4.93 


1.58 


2.83 


0.51 


0.12 


0.31 


9.13 


B04 


1.98 


2.64 


1.09 


1.35 


0.44 


0.17 


0.22 


6.14 


B06 


1.59 


2.15 


0.87 


0.89 


0.43 


0.21 


0.18 


4.94 


BOS 


1.36 


1.58 


0.75 


0.66 


0.42 


0.24 


0.15 


4.40 


BIO 


1.21 


1.08 


0.66 


0.52 


0.42 


0.27 


0.14 


3.71 


B12 


1.09 


0.90 


0.60 


0.43 


0.42 


0.30 


0.13 


3.31 


B14 


1.01 


0.79 


0.55 


0.37 


0.40 


0.33 


0.12 


3.08 


XOl 


6.03 


5.78 


3.32 


2.63 


2.03 


0.27 


0.14 


18.79 


X06 


1.01 


0.89 


0.55 


0.43 


0.34 


0.27 


0.14 


3.07 


XIO 


0.60 


0.57 


0.33 


0.25 


0.20 


0.27 


0.14 


2.78 


X15 


0.40 


0.48 


0.22 


0.16 


0.14 


0.27 


0.14 


1.79 


X20 


0.30 


0.45 


0.17 


0.12 


0.12 


0.27 


0.14 


0.86 



''Computed for post-shock conditions without ambipolar diffusion (see Equations (|52p and (|55|l 'l. 
'^Defined as when the slope of the F vs. time curve drops to 25% of its maximum value. 
'^Lcorc = ^/{n) in "candidate core" region with enhanced n/B. 
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Fig. 9. — Time evolution of the normalized central mass-to-flux ratio F in the shocked gas. The pa- 
rameters are no = 200 cm~^, vq = 5 km/s, Bo = 10 /xG, and Xio = 5 (solid line), with modifications 
as noted in the key. 
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measure the mass-to- flux ratio 'E/{By) at time t = 2fAD! and define this mass-to- flux ratio inside 
the central peak (multiplied by 27r\/G) as Fgnai- 

Table [2] shows the predicted values of tAD and Ffinai from Section 16. 2^ as well as the sim- 
ulation results for these quantities. The measured ambipolar diffusion time scale is ~ 0.3 — 
3 Myr. Our model predicts the ambipolar diffusion time scale very well: the RMS value of 
{tAD, pred " ^AD, sim) /Ud, sim is 0.19, and the range of {tAB, pred - *AD, Sim) Aad, Sim is -0.42 to 
0.28. The measured mass-to-flux ratios deviate from predicted values somewhat more, with a 
range of (Tfinai, pred - Tfinai, sim) f^&nai, sim "0.44 to 0.49, and RMS value 0.28. The typical size 
-^^core = N/{n) of the region with enhanced mass-to-flux ratio at time 2tAD is ~ 0.2 — 0.6 pc (Tabled 
column 6). 

In most of our simulations, the mass-to-flux ratios are higher than 0.6 (see column 5 of Table[2]), 
with some cases (NIO, V12, B02, B04, and XOl) reaching r2t^n > 1. Recall that we assumed the 
effective length of the system is comparable in all directions when we define T in the candidate 
core material (see Section 15. ip . This means that the real mass-to- flux ratio would differ from our 
measured ^2tAD a factor Ly/Lx- Since cores may have axis ratios ~ 2 : 1, the measured ^2tAD 
may be underestimated by a factor up to ~ 2. Therefore the fact that almost all models have ^2tAo 
close to 1 shows that C shock transients may lead to supercritical cores quite frequently. 

The models with T2tj^Q > 1 confirm our prediction that small values of Bq are crucial for 
forming supercritical cores (otherwise uncommonly high neutral density/inflow speed or extremely 
low ionization fraction may become necessary). Given the limits on physical conditions in clumps 
within GMCs (see discussion in Section 16. 2p . prompt supercritical core formation would prefer- 
entially occur if the inflow direction is aligned relatively close to the magnetic field. A study of 
oblique shocks using similar analysis to that in the previous section is performed in Appendix [Aj 
where we show that the transient behavior of C shocks is insensitive to the component of magnetic 
field parallel to inflow velocity, so that our 1-D model is qualitatively applicable in cases with more 
general geometry. 

We also list the value of T2tA^/^BE in Table [2] for all our numerical models. In most cases, 
this ratio is greater than 2, and the average value is ~ 5.4. This means that for essentially all 
reasonable parameters, transient ambipolar diffusion will be important in enhancing the mass-to- 
flux ratio for forming cores. Since r2tA.D close to 1 in many situations, and r2tAD/rBE is large, 
transient ambipolar diffusion during core formation clearly plays an important role in setting the 
stage for subsequent core evolution. 

7. Summary 

Ambipolar diffusion is an important phenomenon in interstellar clouds, which are strongly 
magnetized in the sense ua ~ ^» c^, but are poorly ionized. Supersonic turbulence creates 
shocks, but ambipolar diffusion between ions and neutrals spreads these shocks out. The thickness 
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of C-type shocks depends on the inflow velocity, density, the magnetic field strength, and the 
ionization fraction. Although C shocks are normally studied in the steady-state limit, their early 
transient development is quite interesting. During this transient stage, the central compression 
of neutrals is strongly enhanced because they are effectively "unmagnetized." The time and space 
scales of these transients make them important to the structure and dynamics within GMCs. The 
transient duration is comparable to the drift time across the C shock thickness (~ 0.1 — 1 pc), 
typically ~ 0.1 — 1 Myr for GMC conditions. 

For star formation, ambipolar diffusion is usually analyzed in the context of slow evolution 
leading to gravitational collapse in magnetically-supported clouds. However, our results show that 
since neutrals can stream through field lines in shocks because of ambipolar diffusion, magnetically- 
supercritical cores may form due to C shock transients. During the transient, strong central com- 
pression of neutrals enhances n/B compared to steady-state values. If the compression and duration 
of the transient are sufficient, the central post-shock region may become supercritical and collapse 
gravitationally to make a prestellar core before it re-expands. 

For both the traditional picture of supercritical core formation and the scenario we propose, 
the magnetic field remains relatively stationary while the neutrals move inward, within the high 
density regions. For the traditional picture, the inward neutral motions are due to small-scale 
self-gravity within the core. For shock-induced core formation, the inward motions of neutrals 
owes to large-scale converging supersonic flows within GMCs (which may ultimately be driven by 
large-scale self-gravity within the cloud). 

Transient ambipolar diffusion is particularly important because without it, post-shock regions 
in GMCs typically have very small mass-to-magnetic flux ratios. Thus, the regions with the shortest 
gravitational timescales (at high density, due to shocks) would be prevented from collapsing by 
magnetic fields, which are also enhanced by shocks. Our numerical simulations show a peak in the 
mass-to-flux ratio, produced by transient ambipolar diffusion. For strong shocks (fo/wA,o sufficient) 
and low enough ionization fraction, our results suggest that supercritical cores can be produced. 

Based on our simulation results and analyses, our main conclusions are as follows: 



1. The dominant factors determining the ionization fraction in molecular clouds are ionizing 
cosmic rays and gas-phase recombination. We derive steady-state equations for C shocks 
including ionization and recombination (Equations (j20p and (j29p ). Analyzing the solutions 
of these equations (Fig. [1]) , we find that ionization-recombination equilibrium is generally an 
excellent approximation, and for this regime is mu ch better than the w idely-applied frozen-in 



1/2 i ^ 

condition. For equilibrium ionization, pi oc pn (jNakand 119761 . 119791 ) so that rj 
our notation, and Equation (I27p governs steady C shocks. 



no 
Bo 



1/2 



m 



2. We have solved the steady C shock ODE over a parameter range of upstream neutral density 
10^ — 10^ cm^'^, inflow speed vq = 2 — lb km/s, upstream magnetic field strength 
2 — 15 /iG, and ionization parameter Xio = 1 — 21 (xio is defined in Equations (pSj) — ([26l) ). 
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Using a multilinear fit, we obtain an expression for the C shock thickness (Equation ()45p ). in 
terms of these parameters. We also obtain an analytic expression for the C shock thickness 
(Equation (j43p ). which is in excellent agreement with the numerical result. The dependence 
-^shock oc {voVAfi)^^"^ I {o-Pi,o) Can be understood based on the requirement for momentum 
transfer mediated by ion- neutral collisions. Our r esult for the C shock thickness is comparable 
to previous estimates fe.g. lDraine &: McKedll993ri. although th e parameter dependence differs 
from the case of "frozen-in" ions (IWardlelll99nl : lLi et al.ll2nnd ;i. 



3. During the transient stage of C shocks, the central column density of gas with enhanced 
mass-to- flux ratio initially grows kinematically (Equation (jSSp l. but this slows after a time 
comparable to the ion- neutral drift time ^ad ~ ^shock/adrift across the C shock (Fig. [8] and 
Equation (j60p ). The duration of the transient from our numerical models (see Table [2]) 
is similar to our analytic estimate, 0.1 — 1 Myr for the regime we have studied (see Equa- 
tions (f62]) — (163]) ). Although the present models do not include self- gravity, the duration of the 
transient C shock is comparable to the time needed for prestellar cores to collapse, from both 



obse rvations fe.g. IWard-Thompson et al. 



(e.g. ICong Qstrikeilboilh . 



2007 



Evans et al.ll2009l ) and numerical simulations 



4. Our finding o f rapid initial enhancem ent in density and mass-to-flux ratio is consistent with 
the results of iKudoh &: Basul (120081 ) . for somewhat different parameter regime. Their sim- 
ulations included self-gravity, and they also pointed out that with appropriate parameters, 
collapsing cores may form due to the initial compression. Our work helps to explain the 
physics behind the rapid collapse they identified, and more generally provides insight into 
other numerical studi es of turbulence-accelerated, magnetically-regulated star formation (e.g. 



Li fc Nakamurall2004l ^. 



5. Over the transient time tADi a column ~ 2noWoiAD of "candidate core" material accumulates. 
By taking the ratio with the post-shock magnetic field strength, we can estimate the mass- 
to-flux ratio of this dense material. Equation (167p gives an estimate of the final dimensionless 
mass-to-flux estimate, F, which is similar to numerical measures of Fgnai (Table [2]). The 
relatively hi gh mass-to-fiux ratios we f ind may explain the weak magnetic fields observed in 
dense cores ([Troland &: Crutcheii l2008l ) . Without ambipolar diffusion, the post-shock mass- 

1 /2 

to-flux ratio on the scale of a Bonnor-Ebert sphere (Fbe ~ l-5cs/ (voVAfl) ) would be much 
smaller than the critical value. In contrast, the mass-to-flux ratio in the candidate core 
material produced within transient C shocks is several times larger than Fbe (Table [2] and 
Equation (j72p '). This large enhancement shows the significance of ambipolar diffusion during 
shock-induced core formation. 

6. In transient shocks that can produce Fgnai ^ !> magnetically supercritical cores can form 
and collapse rapidly. Shocks that can reach Fgnai ^ 1 have tAD comparable to the gravita- 
tional free-fall time of the larger-scale cloud (Equation (j69p ). Thus, shock- induced ambipolar 
diffusion is rapid, wherever it occurs. 
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7. The most favorable conditions for forming gravitationally bound cores in cold, turbulent, 
magnetized clouds are strong shocks {vq ^ '^A,o) in regions with low ionization fraction 
{XiO ^ !)• Equation ()A2ip in the Appendix gives the final mass-to- flux ratio for the case 
of oblique shocks; the result is similar to that with the same component of the magnetic 
field parallel to the shock front (Equation ()67p ). Considering realistic conditions in molecular 
clouds, converging flows with Vjnflow -L Bdoud will have relatively low post-shock density and 
post-shock mass-to-flux ratio, and the post-shock gas layers formed will be unfavorable for star 
formation. Cases where Vjnflow and Bdoud are more aligned are more favorable for reaching 
Efinai > 1 (Equation (|7ip). Further observations of the directions of magnetic fields relative 
to observed gas filaments with or without embedded cores will test whether these orientation 
effects are indeed important. If orientation of shocks is in fact important in producing cores 
that can collapse, this may help explain the observed inefficiency of star formation in GMCs. 

While the present models are extremely useful for explaining the phenomenon of transient 
ambipolar diffusion, simulations of more generalized cases are required to support the scenario 
of prompt supercritical core formation in shocks. Three-dimensional simulations of systems with 
oblique shocks, including self-gravity of the gas, would be immediately helpful. In addition, a more 
realistic core-forming environment can be examined by adding nonlinear turbulence to the inflow 
velocity fleld. Further simulations along these lines, together with observations probing density and 
magnetic structure in filaments and cores at different stages, will improve understanding of what 
precipitates star formation. 

We are grateful to the referee for a very thorough and helpful report. This work was supported 
by NASA under grant NNX10AF60G. 

A. Oblique C shocks 

The main text considers a 1-D system with velocities and magnetic fields perpendicular to each 
other, for simplicity. We expect that our results will qualitatively hold for more general geometry. 
Here, we show that under certain conditions, our results can quantitatively be applied to oblique 
C-type shocks. 

In the following, we shall consider a plane-parallel shock in the standard shock frame, using 
the same coordinate system as before. The shock front is in the y-z plane, the upstream flow 
is along the x-direction (vq = vqx), and the upstream magnetic fleld is now in the x-y plane 
(Bcioud = B^flic + Byflf), at an angle 6 to the inflow {Byfi/B^fl = tan 6*). 

For steady, plane-parallel shocks, dt = dy = dz = 0. From the mass and momentum conserva- 
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tion equations for neutrals (Equations (jla p — (jlbp ). we have 

d 



{PnVn,x) = 0, (Ala) 
^ {Pnvl^x + Cs/On) = apiPn (Vj.x " ^'n.x) , (Alb) 
^ iPnVn,xVn,y) = apiPn {Vi^y - Vn,y) ■ (Alc) 



Similarly, the momentum equation for ions and the magnetic induction equation (Equations (|2ap — (pb 
are (with the strong coupling approximation) 



1 d . 



By = apiPn {vn,x - Vi^x) , (A2a) 



Svr dx 
4tt dx 

Vi,xBy - Vi^yBx = const. = voByfl. (A2c 



By = -apiPn {Vn,y - Vi^y) , (A2b) 



Note that B^ = const. = Bxfi in plane-parallel shocks, since V • B = 0. 
By defining our parameters as 



r.x^^. rB^rs, = 4^, (A3) 



PO Vn,x Vi^x Byfi 

and 



M^Mx = 'A ^ = ^ = U^V, (A4) 



Cs ' (3y Snpoc^ 

and applying the ionization-recombination equilibrium (pi cx pn ), Equations (fT2]) and (fT3]) become 

^rl = -C-^M\T (- - -\ (A5) 

dx Vq \rix TnJ 

-(—\+-{r^) = ---r\ (A6) 
dx \ Tn J dx " I5y dx 

Since Equation (jAGj) is the same as Equation (jl3p with /3 — t- (3y, the vs. r„ relation for oblique 
shocks is the same as in shocks with Bx = (Equation (jl6p ). In addition, in the post-shock regime, 
Tnj = fix J. From Equations (|AlcP and ()A2bp . the neutral velocity parallel to the front is given by 

{tb - 1) . (A7) 



1 / VAy,0 



The governing equation in the direction perpendicular to the shock front now becomes 

d -Drl!"^ ( 1 1 



This approaches Equation (j27p only if r^x ~ '^b- 
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It is straightforward to show (|Wardldll99ll ) that 

1 _ 1 tan^ e + (r„rB)"^ + 2 (r^ - 1) {fiyM'^y^ _ 



tan^ 9 + r^^ 



(A9) 



evidently rj^. « for — )- 7r/2. By substituting for tq in terms of using Equation (fT6|l . 
Equation ()A9P gives r^a; in terms of r„. Using this, Equation (jASP may be integrated. Fig. [10] 
shows an example of the C shock structure with varying 9 values. For sufficiently large 0, the shock 
is changed little with respect to the 9 = 7r/2 case. For smaller 9, the structure is quantitatively 
different, but qualitatively similar. 



The final magnetic compression ratio rsj is obtained from Equation ()A9P using rf{9) 



ixj 



rf{9) (tan^e-^) 



tan^ 



2rf(e) 



(AlO) 



From the exact solutions, we know that rf{9) < rj(7r/2) = rj^go, Equation (]A10p therefore suggests 
that for each model, there is a minimum angle between Bq and vq: 



tan 



' 2r/^90 



tan 



1/2- 



(All) 



Since, for a given By^ (or /3y), small 9 corresponds to large upstream magnetic field strength, a 
shock is not possible for very small 9. More practically, this can also be written as 



^^0 



1 — sin 



, cloud" 



sm ( 



or fc' > 



(assuming sin^min ^1/ sin^min) for 



sm I 



V2 



^A, cloud 



(A12) 



(A13) 



where fA,cioud = -Bdoud/ \/4vrpo- The reason for the condition 9 > 9ram is to ensure that the inflow 
is strong enough to produce shocks in the magnetized gas. 



To obtain rf{9), we need to simultaneously solve 



rf{9) 



+ rf{9) + 



BJ 



(A14) 



and Equation (lAlOp . which can only be done numerically. Alternatively, we can also use Equa- 
tion (IA14jl to write (assuming M'^ » r/(0) » 1) 



'^BJ 



^Py[M^-rf{9)] « ^yM (l - ^) 



(A15) 
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Substituting Equation (jAlSP into Equation (jAlOp gives us a quadratic equation for rf{9): 



rfHO) 



fKM^ tan2 e 



+ 



PyM"^ tan2 e 



(3yM 



2M 



Since Ai ^ 1, keeping only Ai ^ terms gives 



(A16) 



rf{d) ^ ^fiyM 



+ 



^KM tan2 B 2M 



+ 1 



(A17) 



This is an analytical approximation for rf{6) (see Fig. llOp. The compression factor rjgg for the case 
with magnetic field parallel to the shock front (tan^ — )■ oo) is ?'/,9o ~ \fWy^ (see Equation (jlSp ). 
Note that for tan0 ^ [2/ (\/^-^)] j '^/(^) ~ '"/.QO- Thus, for all but the smallest angles, oblique 
shocks have similar compression factor to the 90° case with the same -By,o- 

Since Vn/rix is small thorough most of the shock just like rn/rs, we can follow the derivation 
in Section [3.31 to get the formula for the C shock thickness with different rf{9). Equation (jASP for 
the shock structure is 

^)=-Dri'-(,-^). (A18) 

' n / \ 'ix / 



dx 



Tn + 



similar to Equation ()37p . Therefore, as for Equation 
written as ^ 



21), the oblique C shock thickness can be 
4^0 



(A19) 



An example comparing the approximation Equation ()A19P (using rf{9) from Equation ()A17P ) with 
the exact solution is shown in Fig. [TOl 

Regarding the time-dependent behavior of oblique C shocks, we use convergent flow to produce 
shocks in numerical simulations. To see how the component of magnetic field parallel to the infiow 
direction (-B^) can affect the evolution of the candidate core material, we fix the values of ng, vq, 
Xio-, and choose different values of 6 so that = By^ cot 9. Based on our theory, the growth 
rate of column density dN{Yi)/dt is proportional to f inflow (^) = [^f{(^) ~ 1] /^/(^)) which should 
be almost the same for different 9^ since rf[B) ^ 1. The ambipolar diffusion timescale t^D and 
the final mass-to-flux ratio Ffinai) however, should decrease slightly for smaller 9 because of their 
dependence on rf{9). The generalizations of Equations (j6T]) and (f66|) are: 



rfinal(^) 



iAD,90 



/— 2pQVQ X tMi{9) 
ivrvG • ~ i 

rBjBo 



1/2 



final,90 



rf{0) 



1/2 



(A20) 
(A21) 



where we apply rsj ~ rj^go to get the second equation, and Equation ([67|) gives Ffinai.go- Since 
rf{9)/rf gQ is order-unity unless 9 is extremely small, rfinai(^) is close to Ffinai.go in most cases. 
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Fig. 10. — The structure, the final compression ratio, and the shock thickness of an obhque C shock 
with no = 500 cm~'^, By^ = 5 ;uG, vq = 5 km/s, and Xio = 5, as functions of the angle 9 between 
Bq and Vq from 6 = ^min = 16.49° to 6 = 90°. The analytical approximations. Equations ()A16P 
and ()A19p . provide good estimates to the exact solutions. 
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Fig. 11. — The transient behavior and time evolution of the column density and normalized central 
mass-to- flux ratio in the post-shock gas of oblique C shocks with no = 500 cm~^, Byfi = 5 ;uG, 
vq = 5 km/s, and Xio = 5. Though the profile of transient central core differs, the growth rate and 
the final value of F are very similar in each case. 
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The simulation results shown in Fig. [TT] agree with our expectation. The column density grows 
at an identical rate in all cases (though the shape of the central peak differs from one to another), 
and the transition happens slightly earlier for smaller 9. There is no obvious difference between 
the final mass-to-flux ratios in each case, however, since tbj actually decreases for smaller 9 and 
makes FfinaiC^) slightly larger, thus cancels part of the effect from rf{9). 

In conclusion, these tests show that the evolution of C shock transients to make candidate 
prestellar cores is not significantly affected by the component of magnetic field parallel to the 
inflow velocity. 
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